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Abstract: We test two-dimensional TPSA of biphoton light emitted 
via ultrafast spontaneous parametric down-conversion (SPDC) using the 
effect of group-velocity dispersion in optical fibres. Further, we apply this 
technique to demonstrate the engineering of biphoton spectral properties by 
acting on the pump pulse shape. 
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The spectral properties of two-photon light are fully characterized by two-photon spectral 
amplitude (TPSA). It was first introduced for describing two-photon light generated via spon- 
taneous parametric down-conversion (SPDC) from femtosecond pulses and the corresponding 
interference experiments |[Tl|2l. Recently, TPSA has been shown to determine the degree of 



frequency entanglement In particular, the technique of heralded generation of pure single- 
photon quantum states relies on testing TPS A H. 

Experimental measurement of TPSA has been first performed by registering the distribution 
of coincidences between signal and idler photons as a function of frequencies selected in the 
signal and idler channel (joint spectral distribution) |l5]|6l|2l- For biphotons generated from a 
cw pump, when TPSA depends on only a single frequency argument, an alternative technique 
has been suggested for its measurement, based on its spreading in optical fibres |l8]|9l- This 
technique, in particular, enabled observation of various Bell states contained within the width 
of the TPSA JTOl . 

In the case of SPDC pumped by short pulses, TPSA depends on both signal and idler frequen- 
cies, as signal and idler photons are not any more delta-correlated in frequency. Propagation of 
pulsed two-photon light through optical fibre has been studied in Ref. UT], and it was suggested 
to use dispersion spreading for the spectroscopic study of conditionally prepared single-photon 
states. In an earlier paper |[T2l . a similar technique was used for the study of single-photon states 
prepared via cw SPDC. 

In this paper, we show that propagation of two-photon states through a dispersive medium 
(for instance, an optical fibre) provides more general information. Namely, it provides a 
frequency-to-time two-dimensional (2D) Fourier transform of the TPSA and this enables one 
to characterize TPSA by measuring time intervals. Further, we apply this method to a proof-of- 
principle demonstration of the engineering of spectral properties of biphoton entangled states. A 
very efficient technique of measuring two-dimensional coincidence time distribution has been 
suggested in Ref. ifTTI . Here, we use a simpler method, based partly on time selection and partly 
on frequency selection. Despite its simplicity, our method allows us to observe such effects of 
pulsed TPSA as asymmetry with respect to signal-idler exchange and interference structure. 

For SPDC from a short-pulsed pump, the two-photon state has the form |IT1|2| 



where {(Os) , a'' (cOi) are creation operators of the signal and idler photons. The TPSA 0),) 
depends both on the pump spectrum and on the phase matching conditions in the crystal 111. It 
is convenient to introduce the deviations of the signal and idler frequencies 0)^,0),- from exact 
phase matching, which, for simplicity, we will consider as frequency-degenerate at frequency 
COo'. cOs = (Oo + COj = Of) + The two-photon time amplitude (TPTA), whose physical 
meaning is the probability amplitude of the signal photon registered at time fj and the idler 
photon at time f,-, is the 2D Fourier transform of TPSA: 



In a dispersive medium, as, for instance, an optical fibre of length /, each of the photon 
creation operators from Eq. ^ acquires a frequency-dependent phase ||9l that can be attributed 
to the TPSA. As a result, the TPSA changes: 



where A:'.', are second-order derivatives of the dispersion law for signal and idler photons. Trans- 
formation ([3]i, as is well-known in the diffraction theory or in the theory of dispersive spreading 
of short pulses, at a sufficiently large length / ('far-field zone') leads to a Fourier transformation 
from frequency to time. Note that this transformation will occur in both £2j and . In a fibre 
without birefringence, A:',' = k'/ = k". As a result, the TPTA after the fibre has the same shape 
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as TPSA, with the only difference that the frequency arguments are replaced by scaled time 
arguments: 

F(f,,f,) oc F{^,,Q.i)\a=,^ik"LQ.i^tJk"i- (4) 
The probability of the pair to arrive at times f,,fi is given by the squared modulus of (|4|. 
Distribution |F(fi ,f,)p can be measured directly if arrival times of the signal and idler photons 
are registered separately, using, for instance, the pump pulse as a trigger fn\. This is not nec- 
essary for our purposes and we measure only delays between the signal and idler photons, in 
a simpler standard START-STOP technique. This way, the information about the arrival times 
of signal and idler photons separately is discarded. The situation can be explained by passing 
from coordinates fj , f,- to a 45° rotated frame t± = [ti ± fj)/ -\/2. Distribution of the coincidence 
counting rate versus the delay between signal and idler arrival times, f, — f,-, gives the integral 
of the TPTA square modulus w.rt. f+, which corresponds, with the scaling £l± = t±/k"l, to the 
frequency distribution of the form 

^(n_)meas-/dn+|F(i2+,n_)|2, (5) 

where, similarly, £l± = (H; ± £1^) / \/2. 

Distribution ^ provides the projection of onto the axis Q. (Fig. 1). Additional 

information about TPSA can be obtained by measuring the same distribution with narrowband 
filters inserted in front of the signal or idler detectors. For infinitely narrow filters at frequen- 
cies iljo or the resulting distributions will be / di2+|F(i2i.o,ii()|' or / d£2+|F(i2j,n,o)P, 
respectively. As a result, the signal-idler time delay distribution will show the projection on the 
£1 axis of the TPSA intersection with the ii, or il, axis (shown by bold red lines in Fig.l). 
The three distributions taken for the cases of (i) no filters inserted, (ii) a filter inserted into the 
signal channel, and (iii) a filter inserted into the idler channel will provide, in an easy way, basic 
information about the TPSA. For the TPSA shown in Fig. la, these distributions are plotted in 
Fig. Id. From the widths of distributions (ii) and (iii), Ai2, and Aily, one can retrieve the tilt 
a of the TPSA as tana — AHi/An, (Fig. la). From any of these two widths and the width of 
the unfiltered distribution AH, the degree of frequency entanglement 131 can be calculated as 

D _ An 1 ^ AH 1 

A^j 1+cota An, 1+tana' 

In other words, the three types of coincidence distributions (cases (i), (ii), and (ii)) provide 
the projections of the TPSA and of its cross-sections. These measurements also allow one to ob- 
serve interference structure in the shape of TPSA. Such structure can be introduced deliberately, 
to 'engineer' the spectral properties of biphoton light and, in particular, to get entanglement be- 
tween multiple pairs of frequency modes il3|. The easiest way to create a complicated shape 
of TPSA is to split the pump pulse in two and hence to introduce modulation into the pump 
spectrum lfT4ll . Typical shapes of TPSA in this case are shown in Figs 1 b,c. In general, the 
interference structure manifests itself in all three cases (i) (no filters), (ii) (narrowband filter in 
the signal channel), and (iii) (narrowband filter in the idler channel). At the same time, TPSA 
modulation with larger period (Fig. lb) is more likely to be seen in the unfiltered distribution 
than in filtered ones, since the TPSA projection will resolve the structure (Fig.le, blue dashed 
line) while its cross-section will not (Fig.le, black dotted line and red solid line). In the case 
of a small modulation period (Fig.le), unfiltered distribution will have the structure 'smeared' 
(Fig. If, blue dashed line) while filtering with sufficiently small bandwidth will reveal the struc- 
ture (Fig. If, black dotted line and red solid line). 

In our experiment (Fig. 2), type-II SPDC was generated in a 5 mm BBO crystal by means of 
a femtosecond-pulse pump with the wavelength 404 nm and bandwidth 2 nm. The pump was 
focused into the crystal by an / = 75 cm lens. After the crystal, the pump was eliminated using 
a mirror and an RG filter, and the SPDC radiation was coupled into a 500 m of S630 (Thor- 
labs) optical fibre with demagnification about 100, by means of a 20x objective lens. After the 
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Fig. 1. (color online) Squared modulus of TPS A plotted in coordinates Q.± = {Q.j±Q.s)/\/2 
(a,b,c). The TPSA was calculated for the case of a 5 mm BBO crystal pumped by (a) single- 
pulse pump and two pump pulses separated by (b) 520 fs and (c) 1 .75 ps. By registering 
coincidence counting rate as a function of signal-idler delay, after transmitting biphotons 
through an optical fibre, one retrieves the projection of TPSA onto the Q.^ axis (figures 
d,e,f, blue dashed lines). If a narrowband filter is inserted into the signal (idler) channel, 
projected is not the TPSA but its cross-section along the fl, (fl,) axis. The corresponding 
distributions are shown in figures d,e,f (black dotted lines and red solid lines, respectively). 

fibre, the beam was polarization-corrected with the help of two retardation plates and sent to a 
polarizing beamsplitter (PBS). Further, horizontally and vertically polarized photons were reg- 
istered, respectively, by detectors 1 and 2 (MPD single-photon counters, jitter time 45 ± 5 ps). 
Distribution of the delay time between signal and idler photons was measured by means of a 
time-to-amplitude converter (TAC) followed by multi-channel analyzer (MCA). The resulting 
time distributions for the cases of (i) no filtering, (ii) a 1 nm filter inserted into the signal (ordi- 
narily polarized) beam, and (iii) a 1 nm filter inserted into the idler (extraordinarily polarized) 
beam are shown in Fig. 3a (respectively, by blue triangles, black squares, and red empty circles). 
For convenience, all dependencies are normalized to unity. The results clearly demonstrate the 
asymmetry of TPSA with respect to signal and idler frequencies, as the same filter inserted into 
signal and idler channels leads to different time distributions. From the experimental distribu- 
tions, the tilt of the TPSA is 72 ± 0.5°, while from theoretical calculations (Fig. Id), its value is 
73°. For the degree of entanglement R, the experimental distributions give a value of 1 .3, while 
the theoretical distributions lead to R ^ 2.4. The disagreement is mainly caused by the finite 
bandwidth of the filter used (1 nm), which broadens the (ii) and (iii) distributions but does not 
change the (i) one, and hence reduces the R value. 

Lines in Fig. 3a are theoretical distributions calculated for all three cases. Insertion of nar- 
rowband filters was taken into account by multiplying TPSA by frequency-dependent Gaussian 
functions describing the filter transmission band. When calculating time distributions, we as- 
sumed the fibre GVD to be k" = 4.3 • lO^^^s^/ cm. This value was found by transmitting a pulse 
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Fig. 2. (color online) Experimental setup. Second harmonic of a Ti-Sapphire laser, after 
being cleaned from the first harmonic radiation by two prisms, is focused into a BBO 
crystal. The pump spectrum is controlled using a spectrometer (SM). After the crystal, the 
pump radiation is cut off by a mirror (M) and a filter (F), and the down-converted light is 
fed into 500 m of standard optical fibre. After the fibre, polarization is corrected by means 
of a half-wave plate (HWP) and a quarter-wave plate (QWP). A polarizing beamsplitter 
(PBS) separates signal and idler radiation and sends each beam to a single-photon detector. 
Distribution of the photon arrival time intervals is measured by means of TAC and MCA. 
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Fig. 3. (color online) Measured (points) and calculated (lines) distributions of the delay time 
between signal and idler photons, for the cases of no filters inserted in front of detectors 
(blue triangles, blue dashed line), 1 nm filter inserted into the idler channel (red empty 
circles, red solid line) and 1 nm filter inserted into the signal channel (black squares, black 
dotted line). Figure b shows the point-spread function of the method. This function was 
measured by registering time delay distribution without the fibre. 
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Fig. 4. (color online) Spectra of the pump with the modulation created by a BBO crystal 
inserted into the pump beam. Left and right figures correspond to different tilts of the crystal 
leading to different phases of the modulation. 

from the Ti-Sapphire laser first harmonic through the fibre and measuring its autocorrelation 
function width. It is very close the GVD of a standard fused silica fibre 11 5 II . The finite resolu- 
tion of the method was taken into account by registering the coincidence distribution without 
the fibre inserted (Fig. 3b). The measured distribution, with FWHM equal to 90 ps, was used as 
a point-spread function, and the theoretical distributions in Fig. 3a were calculated as its convo- 
lution with the 'raw' spectra, given by Eqs. ( I4l5l l. The theoretical curves for the filtered cases 
are in perfect agreement with the experimental data, although the only fitting parameter is the 
height of the distributions (all are normalized to unity). The 10% difference in the experimental 
and theoretical widths of the unfiltered distribution can be due to the spectrum restriction by 
optical elements and the detectors sensitivity curve. 

In the second part of our experiment, we 'engineer' a TPS A with modulated shape and ob- 
serve this shape using our method. For this purpose, we modify the pump spectrum using a 
birefringent material with the plane of the optic axis oriented at 45° with respect to the pump 
polarization. This method, first suggested in Ref. |[T4l . allows to split the pump pulse in two 
temporally separated pulses with orthogonal polarizations, oriented at ±45° to the initial pump 
polarization. The spectrum of the pump pulse is then modulated by a harmonic function, the 
modulation period scaling as the inverse time interval between the two pulses. By introducing 
modulation into the pump spectrum, we also change the TPS A (as shown in Fig.Ib,c). As a 
birefringent material, we used a BBO crystal of length 1 mm and the optic axis oriented at 45° 
to the pump beam, which led to a 350 fs separation between the two pulses. Note that the crystal 
did not produce SPDC radiation as its orientation did not satisfy the phase matching conditions. 
After projecting the pulses onto the same polarization direction, we obtained the pump spec- 
trum modulated by a sine/cosine function, depending on the phase between the pulses, which 
could be aligned by slightly tilting the crystal. Fig.4 shows the pump spectra with the modula- 
tion phase fixed at (a) and K (b). Further, we only focused on the case of the n phase as the 
shape is more recognizable. 

Experimental and theoretical distributions (points and lines, respectively) are shown in Fig. 5. 
In the case of no filters inserted (blue triangles, blue dashed line) interference structure is clearly 
seen. For filters in the signal or idler channels, we obtained no interference structure, as pre- 
dicted by theory. The plot shows the case where a 1 nm filter is inserted into the signal channel 
(black squares, black dotted line). 

TPSA with the structure shown in Fig. Ic can be obtained by inserting into the pump beam 
a 5-mm BBO crystal with the plane of the optic axis oriented at 45° to the vertical axis. We did 
this in the experiment and indeed, the interference structure in the unfiltered distribution was 
smeared, as one can predict from Figs. lc,f. At the same time, to see the interference structure 
in the filtered distributions one should perform much more narrowband filtering than with the 
interference filters we used. Filtering with 1 nm filters did not reveal the interference structure. 
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Fig. 5. (color online) Measured (points) and calculated (lines) distributions in the case of 
the pump spectrum shown in Fig. 4b. Blue dashed line and blue triangles correspond to the 
unfiltered case, black dotted line and black squares, to the case of a 1 irai filter inserted into 
the signal channel. 

In conclusion, we have performed an experiment aimed to demonstrate a technique for study- 
ing the TPSA of ultrafast (femtosecond-pulsed) two-photon states. By measuring the distribu- 
tion of the time delay between signal and idler photons, we are able to retrieve the projection of 
TPSA on the £2+ axis. With the help of narrowband filters inserted into signal and idler chan- 
nels, one can measure the tilt of TPSA and the degree of entanglement. Further, we used this 
technique for showing the effect on the shape of TPSA of the pump spectrum modulation. This 
last result paves the way, in view of application to quantum technologies! 161, for engineering 
the biphoton TPSA, as, for example, the generation of controlled optical frequency combs |fT3]| . 
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